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Abstract 

^D ' We generalize the single-field consistency relations to capture not only the leading term in the 
squeezed limit — going as 1/g^, where q is the small wavevector — but also the subleading one, going 
as l/(?^. This term, for an (n + l)-point function, is fixed in terms of the variation of the n-point 
function under a special conformal transformation; this parallels the fact that the 1/q^ term is 

'nT , related with the scale dependence of the n-point function. For the squeezed limit of the 3-point 
function, this conformal consistency relation implies that there are no terms going as l/q"^- We 
verify that the squeezed limit of the 4-point function is related to the conformal variation of the 
3-point function both in the case of canonical slow-roll inflation and in models with reduced speed of 
o ' sound. In the second case the conformal consistency conditions capture, at the level of observables, 
1^ ' the relation among operators induced by the non-linear realization of Lorentz invariance in the La- 
grangian. These results mean that, in any single-field model, primordial correlation functions of C 
are endowed with an S'0(4, 1) symmetry, with dilations and special conformal transformations non- 
^j^ I linearly realized by C- We also verify the conformal consistency relations for any n-point function in 

Un ' models with a modulation of the inflaton potential, where the scale dependence is not negligible. Fi- 

^^ , nally, we generalize (some of) the consistency relations involving tensors and soft internal momenta. 

rn 
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O ■ 1 Introduction and main results 



The initial conditions of the Universe are described by the correlation functions of the curvature 

perturbation ( and of primordial tensor modes 7. These correlation functions are expected to be 

^ I approximately scale-invariant, and indeed this matches what is observed for the 2-point function of 

H ■ (^. In many physical systems dilation invariance is promoted to full conformal invariance, so that 

■ ■ ■ it is natural to ask what is the role of conformal transformations for the statistics of primordial 

perturbations. In this paper we study this question focussing on (generic) single-field models. 

First of all, let us understand why dilations and full conformal invariance have a different status 
in inflation. Inflation takes place in an approximate de Sitter space 

a maximally symmetric space whose isometry group is 5*0(4, 1), which coincides with the group of 
conformal transformations in 3-dimensional Euclidean space. The time-evolving inflaton background 
is homogeneous and rotationally invariant, so that translations and rotations are exact symmetries 



of the whole system and show up — in an obvious way — at the level of correlation functions. This 
accounts for 6 out of the 10 generators of S'0(4, 1). The dilation isometry 

t^t- H'^ log X x^Xx (2) 

is not, on general grounds, a symmetry of the inflaton background. Indeed the solution — which we 
can always assume to be of the form = i by a suitable field redefinition — is not invariant under 
a time shift. However it is quite common that there is an additional (approximate) symmetry, the 
shift symmetry (/>—)•(/) + c. In this case the inflaton background is invariant under a combination 
of a dilation and a shift of the flelcu. This unbroken symmetry is the origin of the approximate 
scale-invariance of the inflationary correlation functions [2j . In Fourier space dilations act as 

Ck ^ >^''%x ' (3) 

so that the 2-point function is constrained to be of the form 

{Cj:^Cj:J = {27rfS{h + k2)^F{he-^') . (4) 

Since C becomes time-independent when out of the Hubble radius, the function F must become 
a constant in this limit. Therefore the spectrum must be of the scale-invariant form. The same 
argument implies that the n-point functions of C are product of the momentum-conserving 5- function 
and a homogeneous function of degree 3(1 — n). 

Besides translations, rotations and dilations, de Sitter possesses three additional isometrics 

t^t~2H-^b-x , x'^x'~b\~H-^e-~^^^ + x^) + 2x\b-x), (5) 

parametrized by the inflnitesimal vector b. Notice that these symmetries at large times, t -^ -foo, act 
on the spatial coordinates as special conformal transformations. The reason why these symmetries 
are not so popular in the study of inflation is that, usually, correlation functions are not invariant 
since the inflaton background (p = t breaks these three symmetries. On the other hand, these 
symmetries are explicit when the source of perturbations does not couple (or couples only slightly) 
with the inflaton: this is the case of tensor modes [3j and of scalar perturbations produced by a 
second field [H [S]. One may wonder whether it is possible, like in the case of dilations, to impose 
an additional symmetry on the inflaton Lagrangian in such a way that some combination of the 
isometrics ([5]) and this internal symmetry is preserved. This would require that the inflaton action 
is invariant under 

(j)^ <j) + b-x , (6) 

so that it is possible to cancel the variation induced on cj) by the isometrics (0). This is a Galilean 
transformation [6] with a spatial transformation parameter. It is certainly possible to write an 
inflaton action which is Galilean invariant [3 [8] , but it is well known that this symmetry is only well 
defined in the Minkowski limit and broken by a curved background [9j . This implies that correlation 
functions, which probe the theory at a scale comparable with the curvature, cannot be invariant 
under the isometrics dS]) g. 



^One can also consider the case in which the system is only invariant under a discrete (and not continuous) 
shift symmetry, see [11. 

^In exact de Sitter and in the decoupling limit, interactions of the form (Di^)" are invariant under spatial 
Galilean transformations. However, it is easy to check that all these couplings are redundant, i.e. they can be 
set to zero by a proper field redefinition. This field redefinition is of the schematic form </)—>(/)+ (Dt/))""^ and 
thus vanishes on large scales, giving trivial cosmological correlation functions. Correlation functions are not 
trivial on short scales and are obviously 50(4, 1) invariant because they are produced by a de Sitter invariant 
field redefinition. 



However, this is not the end for special confornial transformations in single-field inflation. Indeed, 
in the case of dilations we know that when we depart from exact shift symmetry (or even exact de 
Sitter geometry) the variation of the n-point function of C under dilation — which is now not zero — is 
related to the squeezed limit of the (n + l)-point function [lOl [TTl |T2]. The purpose of this paper 
is to generalize this argument to the special conformal transformations and to derive conformal 
consistency relations for single-field models. 

We will start in Section [2] to consider the simplest case, when the background metric can be 
approximated with de Sitter space and the inflaton perturbations can be studied in the decoupling 
limit, neglecting the effect on the metric. Moreover, we will assume exact scale invariance, i.e. exact 
shift symmetry on (j). Within these approximations the metric is de Sitter and thus invariant under 
([5]). Special conformal transformations are non-linearly realized, since eq. ([6]) induces a linear depen- 
dence on X of the background solution (similarly, dilations induce a homogeneous shift, (p ^ (p + c, 
of the background). This implies that, in the same way the effect of a homogeneous long- wavelength 
mode on an n-point function is given by a dilation, the effect of a constant gradient is given by a 
special conformal transformation. We thus obtain a relation between the variation of an n-point 
function under a special conformal transformation and the squeezed limit of the (n -|- l)-point func- 
tion, in particular with the term going as l/g^ as the momentum q goes to zero. This generalizes the 
consistency relations for dilations, which relate the variation of an n-point function under dilations 
with the 1/q^ divergence of the (n -|- l)-point function. As the conformal variation of the 2-point 
function trivially vanishes, we recover the known result that the 3-point function does not contain 
terms going as l/g^ in the squeezed limit [13]. The first non-trivial relation is between the 3- and 
the 4-point function and we verify it for models with a reduced speed of sound. In this case these 
consistency relations capture, at the level of observables, the relation among operators induced by 
the non-linear realization of Lorentz invariance in the Lagrangian. 

In Section [3] and Appendix [B] the conformal consistency relations are generalized beyond the 
decoupling limit, i.e. when perturbations of the metric become relevant. Even though in this case the 
metric deviates from de Sitter, and therefore does not possess the 50(4, 1) isometry, the conformal 
consistency conditions can still be derived, using the fact that a special conformal transformation 
of the spatial coordinates induces a profile of C which is linear in x. In other words, with a suitable 
change of coordinates which acts on the spatial ones as a conformal transformation, one can add to 
a given solution a long-wavelength (" mode including the first gradient corrections. This is similar 
in spirit to the treatment of adiabatic modes done in |24j; we study an example of this technique 
(unrelated to inflation) in Appendix[Al In this way the conformal consistency relations can be applied 
to single-field slow-roll inflation and in fact we are able to check the relation between the 3- and 
4-point functions for this class of models. 

In Section[3]we generalize our consistency relations to include deviations from scale-invariance, in 
which the squeezed limit contains both terms going as 1/g^ and terms going as 1/q'^. As an example 
we verify the obtained consistency relations for all the n-point functions in models with resonant 
non-Gaussianity, where the deviation from scale-invariance cannot be neglected. The impatient 
reader can look at the most general consistency relation, eq. ([M]) . In Section [5] and Appendix O 
we generalize some of the consistency relations involving tensor modes and we explain why this is 
only possible in some cases. Finally, in Section [6] we generalize the consistency relations for the case 
in which one of the internal momenta is small. Conclusions and future prospects are discussed in 
Section [71 



2 The limit of exact de Sitter 

In this Section we consider single-field inflation, in the limit in which the background metric can be 
approximated with de Sitter space, i.e. at zeroth order in slow-roll. Furthermore, in most inflationary 
models one can study the inflaton perturbations as decoupled from gravity and thus take the metric 
as unperturbed. In this Section we assume that this is the case, and work in this decoupling limit. 
We also assume exact scale invariance, i.e. exact shift symmetry on <j). When these approximations 
hold, the metric is given by an unperturbed de Sitter metric, equation ([1]), and is invariant under 
the full isometry group 50(4, 1). 

Without loss of generality, through a field redefinition, one can take the background solution for 
the inflaton field to be </>o(i) = t. Let us parametrize the perturbations around this background as 

(t>{t, X) = t + TT{t, X) . (7) 

The background evolution of the inflaton breaks some of the isometrics of de Sitter: dilations and 
special conformal transformations (eqs. ([2]) and dS])). We are interested in studying the conditions 
that the non-linear realization of these isometrics imposes on the n-point correlation functions of 
the field perturbations. 

2.1 Derivation of the consistency relations in exact de Sitter 

As a warm-up let us consider dilations first. It is easy to see that under a dilation parametrized by 
A, eq. ([2]), the perturbations transform as 

TT{x)^it{x) = TT{x) + H-^logX. (8) 

Since the action is here assumed to be invariant under a shift of the field vr — )• vr -|- c, the background 
is invariant under a diagonal composition of a dilation and a shift. Thus, the n-point correlation 
functions are expected to be invariant under this diagonal transformatiorO 



5d{iT^^ . . . iTi:J = {2TTf5{ki + --- + kn) 



n 

3(n - 1) + ^ 4 • 4„ 

a=l 



K^,---^a;J' = 0' (9) 



where (•••)' denotes, here and in the following, a correlation function where the Dirac delta (2'7T)'^d{ki + 
■ ■ ■ + kn) of momentum conservation has been dropped. This implies that the correlation functions 
in Fourier space for large t are homogeneous functions of the momenta of degree — 3(n — 1), that is 



(vr r . . . vr 7* )' = —-. -riTTr . . .TTr )' ■ (10) 

\ aki ak„' Q,3{n-1) \ fci fe„' V y 

The unbroken diagonal symmetry also has another interesting and well-known consequence, 
namely the vanishing of the squeezed limit of a n-point function. In order to see this, take the 
long-wavelength mode vr^ to be approximately space independenlo. Note that equation ([8]) implies 
that we can approximate the effect of such a constant long-wavelength mode on an n-point function 
as a rescaling of the coordinates 

(7r(fi) . . . 1T{Xn))7TL = (^(-^^l) • • • T^i^^n)) , (11) 



•^Note that this is indeed the diagonal transformation; in transforming the n-point correlation function 
we cancel the shift of the field induced by the dilation with an opposite shift. In other words the diagonal 
transformation amounts just to a rescaling of the coordinates in the tt correlation functions. 

'''We are also assuming ttl to be constant in time. This is true out of the horizon in the class of models 
under consideration in this Section up to slow-roll corrections. 



where A = e~^'^^. However, equation ^ tells us that a rescaling of coordinates leaves the n-point 
function invariant. There is then no correlation between an n-point function and a long-wavelength 
mode: the squeezed limit of the (n -|- l)-point function vanishes. 

If we were to drop the assumption of the shift symmetry, the variation of the n-point function 
under a rescaling of the coordinates would be different from zero. The argument above tells us that 
in that case the squeezed limit of the (n + l)-point function would be proportional to the deviation 
of the n-point function from scale invariance (see figured]) 



i-.--k.---^kj"'=HPA<l) 



3(n - 1) + ^ 4 • 4a 



a=l 



i^k. 



TTr 



where Pj^ is the power spectrum of the field perturbations defined as 



(12) 



(13) 



The reader can consult Reference |12] or Section U] for a precise derivation of this result. We thus 
obtain a relation between the squeezed limit of the (n + l)-point function and the dilation of the n- 
point function, given by the non-linearly realized dilation symmetry. We will now see that the same 
happens when considering the special conformal transformations which are spontaneously broken. 
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Figure 1: A constant background mode can he removed by a dilation. A background with a constant 
gradient can be removed by a special conformal transformation. 



Let us start by expanding the long-wavelength mode to linear order in x thuqj 

ttl{x) = ttl{0) + diTrL{0)x\ 



(14) 



^Without loss of generality we are expanding around x = 0. We will see in Section U] that our results do 
not depend on this choice. 



We saw above that the constant piece 7ri(0) can be induced by a rescahng of coordinates which, 
in the presence of a shift symmetry, leaves the correlation functions invariant. We will now focus 
on the linear piece di'TTL{0)x^. Let us write down the transformation of the perturbations under a 
special conformal transformation, eq. ([5]), parametrized by the infinitesimal parameter b 

tt{x) ^n{i) = Tr{x) + 2H-^b-x. (15) 

According to this relation a constant gradient can be induced by a special conformal transformation 
on the coordinates with b = -{l/2)HdTTL{0) B) 

{tt{xi) . . . 1T{Xn))dnL = (^(^l) • • • 7r(x„)) , (16) 

where x* = re* — 6*x^ + 2x*(6 • x). This equation holds in the out-of-the-horizon limit, when one can 
drop the time-dependent piece in the transformation of the x* in the isometry ([5]) and neglect the 
time-dependence of vr. As in the case of the broken scale invariance, the variation of the correlation 
function under such a coordinate transformation is not zero and will therefore induce a relation 
between the variation of the n-point function and the squeezed limit of the (n -|- l)-point function 
(see figure [T|). Let us proceed by computing this variation in Fourier space [3j 

(5b(7r(xi)...7r(x„)) = 



(27r)3---(27r)3^''fci---''fcn^ 




x{-i)y]kl Uxl-2x\{b-Xa) 
a—1 


gj(fei-xiH \-k„-x„) 


•i ^1 d /c„ j(fc^.^j+...4.fc^.^„) 





(27r)3--- (27r)3 

n 

K {-^) ^ b^ [6dl - k^A. + 2K ■ dkjQ {^^ ...hJ- (17) 



a=l 



We can thus compute the (n -|- l)-point correlation function where one of the modes is much longer 
than the other n to linear order in the gradient of the long mode 

n 



a=l 



i^n 



a=l 

where in the second equality we used the fact that a constant mode leaves the correlation functions 
unaffected in the presence of the field shift symmetry. Since within the approximations used in this 
Section one can relate the curvature perturbation in comoving gauge C with the field perturbation 
by C = —Htt, we rewrite this last equation aqj 

iC,<k, ■ ■ ■ CO' '= -2^(1)1' E (e^L - ^^C + 24 • d,^dl) (C,^, . . . CgJ' , (19) 

a=l 



^The long wavelength mode is taken to be time- independent: we are implicitly assuming that its time- 
dependence arises at ©(fc^) and can thus be neglected for our purposes. 

''Here we used the fact that S{q + fci + • • • + fc„) « S{ki + • • • + fc„) up to subleading corrections in q, and 
that the transformation commutes with the Dirac delta (see Section |4|). 
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where P is the power-spectrum of Q. 

In the case of the 2— and 3— point functions, the correlator only depends on the magnitudes of 
the momenta \ka\ = k^, in this case the conformal transformation simplifies to [3] 



^Q-ka 



a=l 



4 d d' 
+ 



,2 



ka Clka t>k„ 



(20) 



In this form it is easy to see that the conformal variation of the 2— point function vanishes (up to 
terms which are singular in Fourier space, see [3] and Section U]), as the two terms in the sum cancel. 
This is obvious from the fact that the 2-point function does not specify any direction with which 
the vector b can contract. This implies that the improved consistency relation just implies that the 
3-point function is suppressed, compared with a local shape, by (at least) 0{q/k)'^ [l3j. Much less 
trivial is the implication of eq. (|19p for the squeezed limit of the 4-point function, as we are now 
going to discuss. 

2.2 Small speed of sound 

In this Section we are going to check the conformal consistency relation between the 3- and 4-point 
functions in inflationary models with reduced speed of sound. These are described by a Lagrangian 
of the form 

S=^ld^x./^{MJ,R + 2P{X,cl))) , (21) 

where X = —■^dfj^^'d'^cl). In order to compute the 3- and 4-point functions we have to expand the 
action up to fourth order, find the interaction Hamiltonian and use the standard in-in formalism. 
All these computations have been done in |151 [16] and here we follow the notation of the first paper. 
Any n-point function can be written as 



(Cfc, • • • CfcJ = {27rfS{h + • • • + k)Pr' n ^-^^"^(^1' •••' ^") (22) 

and the 2-point function is given by 

Kj,Cj; = (2.ym+fe)n^. n = ^^. (23) 



hd 

i=l'^^ 



,3.- -^^^1 ^ 1 ^' 



with e = —H/H and the speed of sound defined as 



"' - P V + 2XPyx ■ ^^^^ 



The amplitude of the 3-point function is: 



^,3, . (^ _ , _ I) M|M , (^ _ ,) ^_ J. g ,.,. , _L g ,.,3 , I ^ ,3 j , ,,,) 

where kt = ki + k2 + k^ and the parameters A and S are related to derivatives of the Lagrangian 
with respect to X: 



A = X'P^xx + ^X'P^xxx 



E = X^P^x + 2X^P^xx . 



(26) 



The 4-point function contains two kind of contributions. One conies from diagrams with a quartic 
interaction and the other from scalar-exchange diagrams. It is easy to reahze that scalar-exchange 
diagrams do not contribute to the conformal consistency relation since they vanish as k\ for A;4 — )• 0. 
This is obvious when a time derivative acts on the soft external leg with momentum k^, as the time 
derivative of the wavefunction gives /c|. The same holds when a spatial derivative acts on the soft 
external leg. Indeed at the relevant vertex we will have external momenta k^ and (say) /C3 and an 
internal momentum —k/i, — k^. Now, depending on which leg the second spatial derivative acts we 
have two diagrams proportional to k^ ■ k^ and — /c4 • {k4 + k^), which cancel each other at leading 
order, leaving terms of order A;| (jfj). This behaviour can be checked in the explicit results of [15tll6j 

We can therefore concentrate on the diagrams with a quartic interaction, whose contribution to 
the amplitude is given by [IS] 



A^il 



Zi _ ^^ ntl kf _ 1 /^3A _ ^ \ fcffc|(fc3 • fc4) / 3(fc3 + fc4) 12fc3fc4 

S S^; kl SU c2+V kf V h + k^ 

1 n _ \ {h ■ k2){k3 ■ h) ( Ei<j k^kj 3kik2k3k4 A I 

32Vc2 V h [^ kl + fc3 Z.^, 



12kik2k3k4 



with fi defined as: 

fi = \x^P,xx + 2X^Pxxx + I 
We are interested in terms linear in k^ which are given by 



+ 23 perm. 

(27) 

(28) 



M 



(4) I 
cont 1^4 -5.0 



4fi-i 



3A 



1 



+ 1 






3^3 

h 



{kl ■ k2){hi -ki 



kik2 + kik^ + /c2^3 , 3kik2k2 



k! 



kf 



+ 2 terms , 



(29) 



where the two additional terms have the same form with ki o k^ and k2 ■^ k^. Notice that since in 
the squeezed limit fci, ^2 and ^3 form a triangle, we can replace (ki • ^2) — >• ^(^3 — ^i — ^i)- Factoring 
out [k^ ■ /C4), we can see that the expression has the same structure as the right-hand side of the 
conformal consistency conditon: 



E(^^' 



4 d 



92 



'■^'^' fc"afc" + M)^^^^l^fc^2^^^3) 



(30) 



Indeed it is straightforward to check that eq. (J19p is satisfied. (To simplify the calculations one can 
choose, without loss of generality, k^ = (1,0, 0). Notice that one has to impose ki + k2 + k^ = at 
the end, after taking all the derivatives.) 

Let us comment on the information carried by the conformal consistency conditions. 



^Notice that it is crucial, for this argument to work, that time derivatives in the propagator are inside 
the T-ordering (or anti T-ordering), so that they simply act on the wavefunctions. In the standard S-matrix 
calculation in Minkowski, one usually has the time derivatives acting outside the 9 functions to restore Lorentz 
invariance, with the contact terms canceling the Hamiltonian terms which are not Lorentz invariant. 



At the level of the Lagrangian, it is well known that operators with a different number of fields 
are related due to the non- linear realization of the Lorentz symmetry [T7] . The conformal con- 
sistency relations contain the same information at the level of observable correlation functions. 
For example, for Cg = 1, the same operator proportional to A (it is the operator {g^^ + 1)^ 
in the language of the effective field theory of inflation [17]), contributes both to the 3-point 
function and to the 4-point function. This is encoded in the conformal consistency relation, as 
the conformal variation of the 3-point function must be matched by a corresponding (squeezed 
limit of the) 4-point function. 

Similarly, if an operator only contributes to the 4-point function and not to the 3-point one, 
then this 4-point function must have a suppressed squeezed limit, i.e. the coefficient of the q~'^ 
divergence must vanish. This is indeed what happens for the operator proportional to /i above 
{{g^^ + 1)^ in the language of the effective field theory of inflation [IZ]). In [18] it was shown 
that it is technically natural to make this operator large and the 4-point function arbitrarily 
large compared to the 3-point function: in any case this happens for any n-point function, its 
q~^ divergence must vanish. 



The conformal consistency relation eq. ()19p states that the q divergence of the 4-point 

1/2 

function is always small, of order P, times the 3-point function. This is not at all obvious at 
the level of the action. Indeed, the operator which reduces the speed of sound gives a 4-point 
function of order P, cj'^ compared with the 3-point function. Moreover, this operator is of 
the form (Vtt)^ so that it naively gives rise to a q~^ squeezed limit. It turns out that, in 
going to the Hamiltonian, the term (Vvr)^ gets a c^ suppression while the leading source of 
the 4-point function (oc c~^) is given by the exchange diagrams which, as discussed above, do 
not contribute to the q~'^ divergence. The same kind of cancellation occurs for the operator 
A7r^(V7r)^ which, in the small Cs limit, would naively give a 4-point which is too large in the 
squeezed limit to satisfy eq. (J19p . Again, this term is cancelled in going to the Hamiltonian and 
gets suppressed by c^, while the leading 4-point function is given by the exchange diagrams. 

Equation (|19p encodes the non-linear realization of the de Sitter isometrics; indeed it looks 
suggestively similar to a Ward identity for non-linearly realized symmetries. Notice, however, 
that cosmologically we only observe correlation functions at late times, so that we unavoidably 
miss some information. For example, the late time power spectrum does not distinguish 
between the de Sitter invariant case and the case with c^ < 1 when de Sitter isometrics are 
non-linearly realized. 



3 Including gravity 



We are now going to see that the conformal consistency relations, which were derived above from the 
non- linear realization of the de Sitter isometry group /SO (4, 1), hold much more generally, basically 
for any single-field model. In this more general case, the group 50(4, 1) arises as the group of 
conformal transformation of 3d Euclidean space. First of all, we are going to show that a long 
wavelength mode can be (locally) created by a suitable change of coordinates, even including the 
first gradient corrections. 



3.1 Adiabatic modes at order k 

We want to show how to construct an adiabatic solution in C-gauge, which is correct not only in 
the homogeneous limit, but also including corrections linear in the gradients. An adiabatic mode is 
locally a coordinate redefinition of the unperturbed solution: we thus parallel what Weinberg [H] 
does in Newtonian gauge and look for the most general transformation of the unperturbed FRW that 
leaves the metric in (^-gauge. As the inflaton is unperturbed in this gauge, the slicing is fixed and we 
cannot redefine the time variable. Since we are not interested in tensor modes, we want our change 
of coordinates to induce only a conformal rescaling of the spatial metric, in this way "exciting" 
only d^. The most general transformation we can do is thus an element of the 3d conformal group 
5'0(4, 1), a different one at each time t. This transformation will turn on (^, that will be given — at 
each t — by a term linear in x and a constant piece: 

C = 2b{t) ■ X + X{t) . (31) 

Now we have to check what are the additional constraints coming from the requirement that the 
solutions we found are the A; — )■ limit of physical solutions [14] . For standard slow-roll inflation the 
linear momentum constraint reads 

dj{H5N - C) = , (32) 

while the Hamiltonian constraint reads 



V2 

(3H^ + H)5N + HdiN' = ^^C + ^HC- (33) 



The coordinate transformations we discussed do not induce 5N (at linear order). Therefore, from 
the first equation, we have to require (" = up to order k'^ if we want our coordinate transformation 
to describe an adiabatic solution at order k. This implies that b and A must be taken as time- 
independent 

C = 2b-x + X. (34) 

The Hamiltonian constraint implies that A^* is of order k. Its spatial dependence would be of order 
k"^ and can thus be neglected within our approximation: A^*(t) just depends on time. We can always 
generate this A^* (t) by a time-dependent translation of the spatial coordinates 

x' = x' + 6x\t) -^ N\t) = -jSx\t) . (35) 

This does not change the metric at constant t, i.e. it does not change Q. The situation is very similar 
for a general model of single-field infiation. In this case the constraint equations are different, but 
from their general form (see for example J13j ) one can still conclude that Q is of order k"^ and A^* of 
order k, which is all that we need in the discussion above. 

Notice that in solving the constraints ()32p and ()33p we have to choose boundary conditions. We 
choose them in a standard way, for example taking C, to vanish at infinity: this does not clash with 
eq. ()34p . as this is only a local approximation of the solution which, once the constraints are solved, 
can be extended to a regular global solution. 

In conclusion, it is always possible to write an adiabatic solution which is valid including 0{k) 
corrections: this is obtained by a (time-independent) special conformal transformation and dilation 
of the spatial coordinates and a time-dependent translation. 
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3.2 Derivation of the conformal consistency relations including 
gravity 

We showed that a long-wavelength mode can be obtained — up to correction of order /c^ — as a coor- 
dinate transformation from the unperturbed FRW. This property holds also in the presence of short 
modes. Indeed, in Appendix [B] we explicitly check, in the case of minimal slow-roll inflation, that 
the cubic action can be obtained with a coordinate transformation of the quadratic action, up to 
corrections which are quadratic in the long-mode momentum. In general, one can constrain the be- 
havior of any {n+ l)-point function in the limit in which one of the momenta q goes to zero. Indeed, 
the effect of this long-wavelength mode on the remaining n-point function can be obtained from the 
n-point function in the absence of the background mode, simply by a coordinate transformation, 
up to corrections of order q^. The coordinate transformation acts on surfaces of constant t as the 
composition of a dilation, a special conformal transformation and a (time-dependent) translation. 
Given that the n-point function is translationally invariant and we are for the moment neglecting the 
(small) scale-dependence, the effect of the coordinate transformation simply reduces to the special 
conformal transformation. 

We can now repeat the arguments that we followed above in the decoupling limit. The infinites- 
imal variation of the n-point function under a conformal transformation with parameter b is given 

by 

Sik.C^k, ■ ■ ■ kj = -^JZ [66-4 - b- 4al + 24 • 4(6 • 4)] {k,k. ■ ■ ■ CO' ■ (36) 

a=l 

This implies, using eq. ()34p . that the formula above with b = V('/2 gives the effect of the long mode 
on the n-point function. To get the (n + l)-point function, we have to multiply by the long (, go to 
Fourier space for it, and finally average over this long mode, exactly as we did above. We get 

(c- c,^,cg, . . . CO' '= -lpi'i)'i'mOk. ■ ■ ■ CO' + o{q/k)' 

" r . . . ^ ^ 1 (37) 

with q'Di = Y^ \^6q- 4 - Q' 4^ + ^4 • dkM' 4)] , 

a=l 

where P{q) is the power spectrum at momentum q. 

This is exactly the same as eq. (|19p that we derived above using the de Sitter isometrics. Now 
we see that it holds in a much more general situation, even when the geometry is not close at all 
to the de Sitter geometry. Notice that in this more general case, the SO {4, 1) group has a different 
interpretation: the conformal group in 3 dimensions. However it is not too surprising that in the 
decoupling limit, when the geometry can be taken as non-dynamical, the equation can be derived 
in a way which has nothing to do with ^, which describes a metric perturbation, but only with 
TT. Indeed in the decoupling limit, the 5*0(4, 1) appears as a symmetry of the action, non-linearly 
realized on vr. 

As discussed in the previous Section, this relation implies the absence of q~^ terms in the squeezed 
limit of the 3-point function. In the next Subsection we use the results above to study the squeezed 
limit of the 4-point function for standard slow-roll inflation. 

3.3 Conformal consistency relation for slow-roll inflation 

Let us verify the conformal consistency relation (j37p for n = 4, i.e. relating the squeezed limit of the 
4-point function with the conformal variation of the 3-point function, for slow-roll inflation with a 
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standard kinetic term. The 3-point function is the classic Maldacena's result [TU] 

H^ 1 



\Qi >'k', ^, 



ki^k2^k3' 



4e2M4n(2A:f) 



(2r? - 3e) ^ fcf + e ^ hk] + £77 E ''^''. 



«7^i 



i>j 



(38) 



where kt = ki + k2 + k^ and 77 is the second slow-roll parameter, related to the inflaton potential V 
hyr] = MlV"/V. 

The 4-point function is the sum of two comparable contributions, one coming from 4-legs interac- 
tions of Q [19] and one from a graviton exchange [20] . The contribution from the graviton exchange 
will not enter in our conformal consistency relation as it decays too quickly in the squeezed limit 
g — 7- 0. This is easy to see from the cubic vertex that enters (twice) in the diagram: 



S 



d'^xdiCdjC'jij ■ 



(39) 



As the graviton is transverse, we can rewrite this vertex moving the two derivatives on the leg with 
momentum q. This implies that the graviton exchange gives a contribution suppressed by two powers 
of q in the squeezed limit compared with a local shape and this does not contribute to eq. (|37p . 
The contact term contribution is given by [19j : 



(40) 



with: 



tj6 1 -,-,-,-, 

P 1 1 \ i J perm 



M4iki,k2,ks,k4^ 



2p u2 u — p + 2A;2 • ^34 + TfTl2CT34 | 



12 "-34 "-* \ "^34 

ki ■ k2 
0-34 I - 
34 V 






1 kl 



kt 



k^k^ 
kf 



/Ci rvn /C4. 



(41) 



W, 



abc 



CTab = ka-kb + kf, , 

Zab = O'abka — O'ba'^b j 
jj^ -, , ka -\- kfj Ikakf, 

ka ~r kf) -\- kc Z[KaKl) + Kb"^c ~r Kakc 



kt 



k! 



+ 



bkakbkc 



(42) 
(43) 

(44) 
(45) 



It is easy to verify that the first term in brackets in eq. (j38|) has vanishing conformal variation 
and thus does not contribute to eq. ([37|) . This must be the case as the 4-point function does not 
depend on r]. The right-hand side of the conformal consistency relation therefore reads 



H^ 



1 



8e^Ml.2q^ll{2kf) V 2 



^Q-ka 



a=l 



4 d 92 
+ 



f^a Oka ok^ 






(46) 



It is straightforward to verify that this indeed reproduces the 1/q^ behavior of the 4-point function 
IQD. 
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4 Including the deviation from scale invariance 

So far we have neglected the deviation from scale invariance, which gives the leading squeezed limit 
contribution ~ l/g^- The inclusion of the tilt is not entirely straightforward, as it is not obvious how 
it will affect the 1/(7^ corrections. Let us go through our arguments once again without assuming 
scale invariance. The starting point is that the effect of a background mode C on an n-point function, 
at zeroth and first order in the gradients of C, just reduces to an appropriate transformation of the 
spatial coordinate: a dilation and a special conformal transformation. In mathematical terms 



(c(fi)...c(xn)y 



C,vc 



(C(^l)---C(^n))o 



with 



Xi 



Xi + C{x+)xi + {Xi - x+)(VC • (Xi - x+)) 



-(x,-x+)M 



(47) 



(48) 



The transformation is done at the midpoint x-|- = (xi + . . . + Xn)/n although, as we will see later, 
this choice does not affect the result. The variation of the n-point function is at linear order in C,: 



5 ^- d \ f d^fci 



dc 
d 



d 



^^'^^rc^'^^ovc 



) J (2vr)3 



d^k, 

W) 



3^^fci---^fc„)o 



(2^) 



3^^fci---^fc„)o 



ifeixiH \-ik„x„ 



ifclXlH \-iknXn 



(49) 



where the tilde on the correlation function indicates the action of the dilation and special conformal 
transformation in Fourier space. To get to the (n + l)-point function, we have to multiply by C{^) 
and average over the long mode. Let us consider the contribution from C, and VC separately. The 
first one gives 



d^q d^ki d^kri 



■P{q) 



-3(n - 1) - ^ 4 • dk. 



(Cg, •••Cj:„)c 



kJO 



(27r)3 (2^)3 ■ ■ ■ (27r)3' 
{2TTf5{ki + ■■■ + k^yk^x^+-+ik„x„e'''^^-^+'^ . (50) 

can be expanded to give transfor- 



(C(x)C(xi)...C(x„))D 



The terms proportional to V^ on the right hand side of eq. 
mations around the origin 

XiiVC • Xi) - -ffVC - XiiVC • x+) - f+(VC • Si) + {x^ • f+)VC + f+(VC • x+) - 24VC • (51) 

The last two terms describe a translation: this does not affect the n-point function which is trans- 
lationally invariant. Analogously, the third-to-last term together with the fourth-to-last describe 
a rotation, which again does not change the correlation function. The third term describes an 
additional dilation and gives a contribution to the (n + l)-point function of the form 



(C(x)C(xi)...C(x„))D 



d^q d^ki d^kr, 



(2vr)3 (2vr)3 " " " (2vr)3 ^^^^ ["^^^ " '^ " ^ ^^'^ ' ^'" 
{2TT)H{ki + ■■■ + 4)e*^^^i+-+^^"^"e^«^(^-^+)(i^- x+) . 



(Cfc, •••Ci:„)c 



kjo 



(52) 



This expression cancels, up to terms of order q'^, the dependence on x^ in eq. (|50p . This shows that 
the choice of the point x_|_ is immaterial. 
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We now come to the action of the special confornial transformation, i.e. the first two terms of 
eq. ()5ip . To derive the action of these in Fourier space we have to convert powers of Xi into derivatives 
with respect to momenta and integrate these by parts. In doing this, some derivatives act on the 
5- function of momentum conservation, as explained in ^, and give a term with a distribution V5 
multiplied by the dilation transformation of the correlation function. This gives 



(C(x)C(xi)...C(x„))D 



d^g d^fci Skn 



(27r)3 (27r)3 ' ' ' (27r) 



P{q) 



-3(n - 1) 



Y^ ka ■ <9fc„ 



(C.A---C, 



kn'O 



q- {2TTfV6{ki + ■■■ + K)e 



ik-ix\~\ \-iknXn iq{x~x+) 



(53) 



Notice that the term q ■ x+ can be dropped as it gives contributions of order q^. If we sum the 
expression above with eq. ()50p we end up with the full momentum conservation (27r)^(5(g + ki + 
. . . + kn)- Besides these singular terms we also have the standard special confornial transformation 
in Fourier space. In total we get to 



(54) 



Notice that the procedure above implies that the scale transformation of the n-point function is 
evaluated exactly with the momenta ki . . .kn that appear on the left-hand side, even if at first order 
in q these will not form a closed polygon. This is irrelevant for the conformal transformation which 
is already at order q. Equation ([M|) represents the most general result of this paper. 







r- ^ 




{(,<l, ■ 


■■4.)"=°-^w 


3(n - 1) + ^ 4 . 4^ + ^qWi 

a 


iCk.---CO' + <^(i/kf 




n 

with q'Di = ^ [eg • 4, - ?• 45|„ + 2ka • dkM- 4 J 

a=l 





4.1 Resonant non-Gaussianities 

Models with a periodic modulation of the inflaton potential [2T] have recently attracted attention, 
as motivated by explicit string constructions of Monodromy Inflation ^22j. For our purposes, these 
models represent an interesting example in which the n-point functions deviate sizably from scale- 
invariance — so that the results of this Section come into play — and in which many n-point functions, 
and not only the bi- and trispectra, can be analytically calculated [231 [1]. 
The (n + l)-point function readqj i23i JJ 

where /„+i is given by the following in-in integral: 

I„.^i = -2Im / ^y(-+i)(<^(^))(l _ iqr,){l - ikir,) ... (1 - ikn7l)e"'''' • 



{Q<k.---CkJ = i^^f^(<i + ki + 



kn 



'n+1 



?]- 



(55) 



(56) 



^This result holds in the decoupling limit and neglecting diagrams with more than one vertex. Corrections 
to the decoupling limit are suppressed by e and they vanish in the limit e — > keeping constant the power 
spectrum of C [1] . Terms with more than one vertex contain extra powers of the amplitude of the potential 
modulation. These corrections eventually become sizeable for n-point functions with large n [23]. However, 
the consistency relations must be satisfied independently by the single-vertex diagrams, as they are the only 
ones proportional to the first power of the amplitude of the modulation. The integral over rj can be explicitly 
done in the saddle-point approximation, in the limit of many oscillations per Hubble time [231 [T]. Since 
the consistency relations must hold independently of this limit, and not only for oscillatory but for general 
potentials, here we prefer to keep the integral form. 
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Here kt = q + Yli ^i ^-i^d y("+^) is the (n + l)th derivative of the potential, evaluated on the 
unperturbed history (j){r]). In resonant models, one starts from a trigonometric function and therefore 
these derivatives are simply sines or cosines. Here we keep the discussion more general, without 
specifying the form of the potential. 

Let us verify our general consistency relation eq. (j54p . First, notice that the action of the operator 
q^Di on the n-point function vanishes. Given that the function only depends on the moduli of the 
wavevectors, the differential operator takes the simplified form (|20p 



^Q-ka 



a=l 



4 d 



+ 



dki 



\^fci ■■■^kj 



(57) 



To see that this expression vanishes, it is enough to compute the action of the operator in brackets 
for one particular momentum ka- The result does not depend on a (indeed the form of the amplitude 
remains the same with an additional —rf that appears in the integral). This means that the sum 
will be proportional to q ■ X^^^^ ka-, which is zero up to q^ corrections. 

This result seems to imply that correlation functions in these models are invariant under special 
conformal transformations. This cannot be true because correlation functions are not scale- invariant, 
due to the features in the potential, and there is no such thing as invariance under special conformal 
transformations without scale-invariance (indeed the commutator of a special conformal transforma- 
tion with a translation contains a dilation). As we discussed above, there is an additional singular 
piece of the special conformal transformation, proportional to the derivative of the delta function 
and to the scale variation of the n-point function — which is not zero in this case — and this was used 
to derive eq. ()54p . The vanishing of the conformal operator in Fourier space tells us that in these 
models, where there is no operator with spatial derivatives that enters at order q in the consistency 
relation, the conformal transformation of the n-point function is just fixed by its dilation properties. 

Let us therefore check that the dilation transformation captures the squeezed limit of the n-point 
function up to corrections 0{q^). In the squeezed limit g — t- 0, the expression for the (n + l)-point 
function gives 



^^T'ki ■ ■ -^kJq^O 



H\ H^f /7\" F2"-4 
2g3 



^y("+i) J|(l - ikir])e'^^'^ , (58) 

-oo V 



up to 0{q'^) corrections. Now we can convert the (n + l)th derivative of the potential using 



-^yW(0(^)) = _S 






(59) 



and integrate by parts the derivative with respect to rj. This derivative can be traded for derivatives 
with respect to the ki to give 



T ^ 3(n-l) + j:^.— 





2g3 


PiQ) 


3(n 



d 



J —oo V 



UU ^kf 



i=l 



3(n - 1) + ^ fci 



dki 



\^fci • • • ^fc„ / ' 



(60) 



as we wanted to show. 
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5 Consistency relations for tensor modes 

In the previous Sections we discussed the generahzation of the consistency relations involving scalar 
modes. There are similar relations involving tensor modes at zeroth order in the gradients |1U] . 
and it is thus natural to wonder whether they can be generalized as well. We will be able to do so 
only in the case of one long-wavelength tensor mode, with the other legs being scalar modes. In all 
other cases, one is not able — at first order in gradients — to create the long mode with a coordinate 
transformation. For example, if one starts from the flat Euclidean metric, a special conformal 
transformation induces a conformal factor with a linear dependence on x, but the same will not 
work if one starts in the presence of a tensor mode. 

Let us consider the (n + l)-point function with one background tensor mode (7^^ ' ' ' Ck ) ^tj)' 
We work in the usual gauge, where the spatial line element is given by a^e^^ {^^)ij dxidxj and 7ij 
is transverse and traceless |10l [3] . We want to induce a gravitational wave in the presence of scalar 
modes only, i.e. starting with a line element a'^e'^''6ijdxidxj. To induce a long-wavelength tensor 
mode, we consider a quadratic transformation of coordinates 

Xi r Xi -p -^ij'^j ~r -Dij^XjX^ , V^-^/ 

where Bij^ is symmetric in the last two indices Bijk = Bi^j (notice this has nothing to do with a 
conformal transformation). With this transformation the spatial line element becomes: 

a'^e'^'^dijdxidxj — )• a^e^^ I 6ij + 2Aij + 4^^ — - — ^Xfc j dxidxj . (62) 

We can identify the terms in brackets with the gravitational wave; imposing the symmetry properties 

of Bijk we find: 

1 1 

^y = 2^^J ' ^*i'= ^ ji^klij - diljk + dj-fik) ■ (63) 

So far we did not put any constraints on jij. From the expressions for Aij and Bijk we see that 
we can consistently impose that jij is transverse and traceless. Now that we have the coordinate 
transformation that induces the long-wavelength tensor mode, we can follow the same logic as with 
scalar modes: 

(C(xi)...C(fn))7 = (C(^l)--.C(^n)), (64) 

where transformation of coordinates is now given by: 

Xi ^ Xi -\- Aij[X^)Xj + IJijk(Xj^)[Xj Xj )[Xk Xf^ ) 

= Xi + {Aij{x+) - 2Bijk{x+)x^)xj + Bijk{x+)xjXk + translation . (65) 

We end up, leaving the details of the calculation to Appendix [Cl with 

i^Kk. ■ ■ ■ CajJUo = -lp.i<i) E ^U'l^kaA^. (Ck, ■ ■ ■ CkJ' 



2 

a 



1„ . .>r- 



4 



-P^{q) }_^ 4^{q) (2kai{q- dkj - {q- ka)dkj dk^^ (C^^ • • • CkJ ■ (66) 



a 



This is the generalized consistency relation for the {n + l)-point function with n scalar modes and 
a long background gravitational wave. The graviton power spectrum is given by P'y{q) = -tpt ■ s- 



''We are using the convention of [10, for the tensor modes. 
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As before, the sum of the momenta entering in the n-point function is not exactly zero and this will 
produce the corrections linear in q coming from dilations. 

We can check this equation in the simple case of the 3-point function [lUj 

In the limit g — )• the amplitude becomes: 

This coincides with the right-hand side of the consistency relation, eq. (|66p . Notice that the second 
line of eq. ()66p vanishes and that linear corrections come from evaluating the 2-point function with 
the sum of the two momenta not exactly equal to zero. 

Let us come back to the question of whether it is possible to generalize the consistency relations 
to cases where we have some short gravitational waves. In this case we start from a line element 
of the form a^e^^ {^'^)ij dxjdxj. If one wants to consider the effect of a long-wavelength mode of C, 
it is still possible to induce it, at zeroth order in gradients, by a dilation of the spatial coordinates. 
But a special conformal transformation induces a conformal factor only if one starts from a metric 
proportional to 5ij and not in general. It is thus not possible to induce a constant gradient of C, 
to generalize the consistency relations. The same arguments apply when considering a long tensor 
mode. The transformation eq. ()63p which induces a long tensor mode including the first gradient 
corrections does not work in the presence of short tensor modes, as now also the part which is 
antisymmetric in ij, which drops out starting from 5ij, gives a contribution. 

All we said can be seen at the level of the action as well. For example, looking at the cubic 
action for one tensor mode and two scalar modes in a slow-roll model 

S = d^x ^jpajijdiCdjC , (69) 

we can see that the effect of 7 is equivalent, including first gradient corrections, to the spatial change 
of coordinates discussed above. A similar analysis does not work for the vertices Cl77 ^^^ 7l77 
when gradient corrections are included. 

6 Soft internal lines 

So far we have considered squeezed limits, in which one of the external momenta goes to zero. 
However, using similar arguments, one can discuss the case when the sum of m (out of n) of the 
external momenta goes to zero. In this case the n-point function is dominated by the exchange 
of a soft — small momentum — intermediate state. The result factorizes in the product of the power 
spectrum of the soft internal line and the m— and (n — m)-point functions, both calculated in the 
presence of the long mode |201 123] , Indeed, the soft mode — scalar or tensor — freezes much before 
the others and behaves as a classical background. Following the arguments of this paper, it is 
straightforward to take into account a constant gradient of the long mode, i.e. the first correction in 
its momentum. In mathematical terms 

+ p,{^) Ei<.<k ■ ■ ■ ^OUiilk^,, ■ ■ ■ ^OU , (70) 

s 
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where 1 < m < n, q = ki + ■ ■ ■ + km and the * on the squeezed-hmit correlation functions means 
that both the delta function of momentum conservation and the long-mode power spectrum P{q) 
are dropped. Notice that, as a consequence of momentum conservation, the vector ^has opposite 
sign in the two terms of the product. We can then evaluate the equation above using our previous 
results eq. ([M|) and eq. ([66|) . 

The simplest example in which an internal line becomes soft is the 4-point function in the limit 
in which two of the external momenta become equal and opposite (counter-collinear limit) |201 123) . 
In the case of scalar exchange we know that the 3-point function {CqCu Cu ) does not have any linear 
corrections. This implies that, in the limit of exact scale-invar iance, we have a q'^ suppression at 
each of the two vertices (this is indeed verified by the explicit results [151 [T6] ) that completely kills 
the g — )• divergence from the power spectrum, so that the result is not dominated by the soft scalar 
exchange. The counter-collinear limit is therefore regular 

i(krCk,Ck,CkJ'g-,0 -- const. (71) 

More interesting is the case of graviton exchange, as in general we have linear corrections. For 
example, plugging in eq. ([TOl) the result eq. ([68|) of the previous Section for slow-roll inflation, we 
have the graviton exchange in the counter-collinear limit {q = ki + k2, ki ~ —k2 and k^ ~ — ^4): 



/A A ^ r \'GE _H^ H^ H^ \^,^*(^,^ t^ ^ ^^'^^'J ^3/^3m /^, , 5 fci • g 5 fcg • g ^ ,^^. 




g3 4efcf 4eA:| Y 4 fcf fef 

We can see that our results allow to capture also the 0{q) corrections in the counter-collinear limit, 
and this is exactly the result that one finds taking the appropriate limit of the full 4-point function 
induced by the graviton exchange [20]. 

We conclude by noticing that as the squeezed limits are important for observables which involve 
a large separation of scales (the scale-dependent bias for example), the same occurs in the case of 
soft internal lines. For example, the recently discussed power spectrum of the CMB /x distortion [21] 
is sensitive to the 4-point function in the counter-collinear limit, and therefore very suppressed in 
all single- field models. 

7 Conclusions and outlook 

In this paper we extended the single-field consistency relations in the limit in which one of the 
momenta — external or internal — goes to zero, to include the first correction in the soft momentum. 
We proved that this correction is related to the variation of the correlator under a special conformal 
transformation. 

These "soft" limits are phenomenologically relevant, as some of the proposed experimental tests 
of non-Gaussianity involve a separation of scales, and therefore probe the primordial correlation 
functions in one of these limits, where our results apply. For example, we showed that in any 
single-field model the 4-point function generally goes as 1/g^ in the squeezed limit (assuming scale- 
invariance), but that the coefficient of this divergence is suppressed by ~ 10~^ compared to the 
3-point function and thus very small. Therefore the scale-dependent bias induced by the trispectrum 
(see [25j and refs. therein) will be suppressed in any single-field model. 

From the theoretical point view our results, besides being a useful check of the calculations, 
fully explore the symmetry of primordial perturbations generated by single-field inflation. In a 
completely model-independent way, in any single-field model primordial correlation functions of C, 
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are endowed with an 50(4, 1) symmetry, with dilations and special conformal transformations non- 
linearly realized by (. Notice that this statement is not as trivial as the simple equivalence between 
a conformal transformation of the spatial coordinates and a change of the scale factor induced by 
C- Indeed this equivalence must be valid during the whole period when perturbations are generated, 
in order to have the S0{4, 1) symmetry at the end of inflation. This holds, of course, only in single 
field models. It is worth noticing that no further extension of these consistency relations for scalar 
modes is possible: at second order in gradients, the long wavelength mode induces curvature and 
therefore cannot be removed by a clever change of coordinates. Similarly to what happens with the 
1/q^ consistency relations, our eq. ([M]) encodes the fact that a long mode cannot be observed locally, 
until one is sensitive to its curvature. 

It is worthwhile to stress that in this paper we switched between two different interpretations of 
the group 50(4, 1). In the decoupling limit it is the de Sitter isometry group while, later on, it is 
the conformal group of 3d Euclidean space. The first interpretation requires exact de Sitter, but it 
can be applied to test fields [S], models with more than one field and to higher spins [3]. The second 
does hold in a general FRW, but it requires a single field model. 

Let us conclude with some comments on possible future directions. Even if we phrased the 
paper using inflation as the key example, our results are more general and can be extended to other 
single-field models like [261 EH ESI [29]. It should also be possible to generalize some of the above 
results to multifield models, at least in the case of exact de Sitter geometry. The inflaton background 
spontaneously breaks this symmetry, so that the variation of a correlation function under the de 
Sitter isometry group should always be connected with the soft emission of an inflaton perturbation. 
However, it is not clear whether this can be simply connected with observable quantities. Similarly 
one should be able to extend these relations to models where the 50(4, 1) symmetry is not the 
isometry group of de Sitter, but a subgroup of the 4d conformal group 50(4,2), with additional 
constraints induced by the non-linear realization of 50(4, 2) [351 [311 [32] ■ 

Work is in progress by other groups to understand the consistency relations in terms of sponta- 
neously broken symmetries [33] , either in terms of Ward identities [34] or using path integrals and 
OPE arguments |35] . 
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A Non-linear metric in matter dominance 

In this Appendix we want to verify in an explicit example that the second-order metric in (^-gauge, 
with one of the modes much longer than the Hubble radius can be obtained simply as a coordinate 
transformation from the flrst-order metric, and that the result is accurate up to corrections of order 
(ki/aH)'^, where k^ is the long- wavelength momentum. We are going to study the second-order 
metric in matter dominance, written in terms of the initial perturbation (^q, assuming no initial 
tensor mode. This has been studied in (^-gauge in [36] . 
The form of the line element is given by: 

ds^ = -N^dt^ + hj{dx' + N'dt)idx^ + N^dt) , hij = a^{t)e^<(6^j + 7^) , (73) 
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N = 


1 




^ = 


-^Co + ^9-2 [(9,Co)' - 3d-^d,dk{d,CodkCo)] 




+ 25a2F3^ 


^{d'Cof + d^Codid'Co + ^{did,Co? 


^r = 


-^d-' [d^Cod'Co - d-^didk{dkCod'Co)] , 


c = 


Co r;.2m^ 


^d,d,id,Cod,Co) ■ 



with 7 transverse and traceless. Setting Ni = diip + Nf we have 

(74) 



(75) 

(76) 

(77) 

Tensor modes 7 will also be induced by the initial scalar perturbations; their explicit form can be 
found in [36], but we will not need it here. We are interested in this metric in the limit in which 
one of the modes in the initial conditions, let us call it Ci, is (for all times) much longer than the 
Hubble radius, keeping only terms which are at most linear in the gradients of this long mode. In 
particular we are interested in the quadratic terms with a long and a short mode, which we call Qs- 
In the limit we are considering these second order terms reduce to 

N, = d,^ + Nf = —±-^djCLdidjCs - ^Cs^^CL , (78) 

C = -^^^^Cl5.C5. (79) 

The linear metric 

ds'^ = -dt^ -diCidtdx' + 0^(1 + 2QL)dx^ (80) 

bH 

can be obtained starting from the unperturbed FRW with a coordinate transformation. One has to 
perform a conformal transformation of the spatial coordinates, which multiplies the spatial metric 
by a constant term and a term linear in the coordinates: in this way one reproduces the spatial 
dependence of C,l up to corrections involving the second derivatives of Ql- To give the correct 
space-time component one has to make a time-dependent translation 

^ = ^ + l^diCL . (81) 

Now we want to apply the same coordinate transformation on the metric at first order perturbed 
by the short mode 

ds'^ = -dt^ diCsdtdx' + 0^(1 + 2Cs)dx^ (82) 

bH 

and reproduce eq.s ()78p and ()79p . Notice first of all that this change of coordinates does not change 
N , which indeed is unperturbed even before taking one of the modes to be super-Hubble, eq. (j74p . 
The first term on the right-hand side of ()78p comes from the time-dependent translation ()8ip of diQs 
in eq. (j82p and the same happens for eq. (j79p . The second term in ()78p comes, on the other hand, 
from the spatial metric, when dx is converted into dt by the time-dependent translation. One can 
also check, from the explicit expression in [36], that tensors are not produced in the limit in which 
we only keep the first gradients of the long mode. 

Notice that we have an additional set of second order terms coming from the fact that the short 
mode Qs is evaluated at the conformally transformed point and this transformation depends on C,l. 
This is not a genuine non-linearity but a choice of initial conditions. 
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B Consistency relations at the level of the action 

In this Appendix we explicitly check that the cubic action for C, in the limit in which one of the 
three modes is much longer than the other two, can be obtained with a coordinate transformation of 
the quadratic action, up to corrections which are quadratic in the long-wavelength momentum. This 
property of the action implies the conformal consistency relations for (^. We will focus on minimal 
slow-roll inflation, whose quadratic and cubic actions are given by [lOj 



5(2) 



d^xdt ea^ 



dC 



(83) 



5(3) 



d^xdt 



^ 1 + 



+ 



a^e^< I i {did,i;did,i; - (9^)') (l " |^) " '^d.i^diCd^^; 



(84) 



where "0 = — tjj + X) d'^X = ^C s-nd we set Mp = 1. Let us consider the coordinate transformation 
that gives rise to the correct linear-order solution of the metric for the long-wavelength mode Cl 



b'x^ + 2x\b-x) + 6x'{t) 



with 



b, 



■iaa |fa'{<) 



2b' 
1^ 



+ 0{e). 



(85) 



(86) 



The slow-roll correction in fe* is necessary to reproduce the x contribution in -0, but we will not 
need it below. Notice that we are here considering a change of coordinates which induces a long 
wavelength mode, while in the main text we consider a change of coordinates which removes the 
background C,l'- this is the source of a relative minus sign. The change of coordinates induces the 
following variation of the quadratic action 



55(2) 



/ 



d'^fdt 



-Qa^e{b ■ x)C^ + 2at{b ■ x){dCf - Aae-^bidiC 



(87) 



Let us compare this with the cubic action at first order in gradients for the long mode. We separate 
long and short modes by writing C, = Cl + Cs and note that C,l = 0{k\). Without loss of generality 
we are expanding C^l around the origin (the calculation for constant C,l and dilations can be found in 
|13j). The relevant terms from the expansion of the cubic action are those linear in C,^. For example, 
the second term in S^^> gives 



6aV^C^fl-4W3ea3aCl. 



H 



(88) 



For the other contributions we have to do some integrations by parts. The first term in S^^ gives 



c 



aeHl + ^ 



2d\ - {dCf) ^ ae^^+fs (^1 + ^) ^_2d\s - {dQsf - 2^,CL^^Cs) 

- aidQsf + 2cMd\s + 2d,QMd,Cs) , (89) 



H 



H 
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where we have done one integration by parts: aCiCsd'^Cs = —a-CLidCs) — at^iCiCsf^jCi- Similarly, we 
can transform the second term in brackets in the following way 

^ -a( - 2d^CL^d,Cs + CUdCsf + eCiidCs?) , (90) 

so that the final contribution from the first term in the cubic action is 

c 



«e^(l + jj) (-25'C - {dCf) ^ -aeCUdCsf. (91) 

Finally, we can focus on the last term in 5*-^'. It is easy to see that the relevant contribution will be 
given by 

^a^CL {didj^Psdidji^s - {d^i's?) - 2a''d^^sd^QLd''il>s - 2a^di^PLdiCsd^iPs ■ (92) 

We can see that the first term has the same form as the second one, once we integrate by parts. 
Moreover, both of these terms are total derivatives. This can be seen by applying the following 
identity 

diid^fdjf) = d^fd.f + difdidjf = d^fd.f + Id.idff . (93) 

We are left only with the last contribution in eq. (|92|) . We are forced to replace d'^ips with d'^xs = ^Csi 
otherwise we again get a total derivative term. At leading order in slow-roll we have diipi = 
— ^diC,L- Collecting all the relevant contributions we get 

5^(3) = j d^xdi Uea^Ul - aeCUdCsf + ^Cs^^Cs^^CL^ ■ (94) 

This coincides with eq. ([87|) using bi = — ^SjCl, as we wanted to show. 

C Derivation of the consistency relation for tensor modes 

In this Appendix we are going to give some details on the derivation of the consistency relation 
for short scalars and a long gravitational wave. We want to compute the average over the long 
gravitational wave ")mn{x) of the variation of the scalar n-point function induced by the change of 
coordinates given in eq. ()65|) . In order to simplify the notation let us define the following: 



Ci0+) = A,j{x+) - 2B,,kix+)x+, P = Y,ki, (Cfc, . . . CfcJ = i27rf5iP)M . (95) 

i 

Under this change of coordinates, the transformation of the n-point function is given by 

n 

5{C{xi) . . . C(^n)) = ^ {Cij{x+)Xajdai + Bijk{x+)XajXakdai) (C(^l) • • • Ci^n)) ■ (96) 

a=l 

We can compute it term by term. The contribution from dilations is given by 

n 
6{Cixi) . . . CiXn)) D ^ Cij{x+)Xajdai{C{xi) ■ ■ ■ C{Xn)) 
a=l 
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■-\-iKYi-Xn 



Here we used the fact that the tensor Cij is traceless and that a derivative on the delta function 
can be rewritten as dp.5{P) leading to a term proportional to Pi5{P). Let us now use the explicit 
expression for Cij-. 

Cij{x+) = -(1 - x+a+)7ij(x+) + - (x+5+7jfc(x+) - x+a+7ifc(f+)) (98) 

and focus on the first term. Averaging over the long background mode ^mnix) we get 

(7..n(x)C(fl)...C(Xn)) d/^|^...|J^(2^)^5(P) (^-E^A.-^) 
s 

Notice that (1 — x'l^d^)e~^'^'^+ = 1 + 0{q'^) which proves that there is no dependence of the result 
on a choice of the point x+ at first order in the small momentum q. The second term of eq. (j98p 
has an additional dependence on x+ but it gives no contribution. After the averaging over the long 
background mode 'jmnix) the result is 

ilmnixKix,) . . . C(X„,)) :^ / (0 1^ • • • §^i^-fm [- E kA,M^ 

X E ''mn{q)^l ikjeUq) - iqie^M) ^^i^^-^^^-^^^-^^S'^-n-^^ . (lOO) 



We can rewrite x'^ in front of the exponent like idq^ and do a partial integration. In the first 
line in the integral nothing depends on q. The action on the second line will produce three terms 
respectively proportional to 

~ E ^^<ik^'mn%) -3 {%4k% - Qi^'jkio)) , 

s 

~ E^-"(^"^§ ii^'ikio) - m]M) , (101) 

s 

~ E ^mniQ)-3 i^jio) + Qjdqketkio) - (^jiio) - Qidq^^jkiq)) ■ 

s 

The first term will vanish when we multiply both sides of the final equation by ^mni^ ^^ order to 
express everything in terms of the amplitudes. It is easy to see that emni'l)(^qk^mniQ) — 0- The 
second term is trivially zero because polarization tensors are transverse. In order to prove that the 
last one is zero as well, let us first note that one can write the polarization tensor in the following 
form 

et;ika) = -^{ef;{ka)+isaefj{ka)) , (102) 
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where 

efj (ka) = ZiZj - uy^, efj {K) = Ziu'^ + Zjuf , (103) 

and unit vectors z and Ua are orthogonal to ka and each other [3] . Using these expressions one can 
derive the foUowing formula for the derivative of the polarization tensor 

dk.etA^) = -^ (heUk) + krefJk)) , (104) 



which one can use to complete the proof. 

We can focus now on the second part of the coordinate transformation. The variation of the 
n-point function is 

n 
6{C{xi) . . . CiXn)) 3 ^ Bijk{x+)XajXakOai{C{xi) . . . C(^n)) 
a=l 

= / ^ • • • ^^^^)' 5;(-^)A,.(f+)<^(P)^^fe«5,„^.a,„,e*•-+■■■+^^^"■- 

= - / n ^(2vr)3 Y, ^B.M^+)ka^ {2d,^AP) + md^) Q^^^Me^k.^^^-^^^^^" , (105) 



where we have used similar tricks as before to get the final result. It is convenient to express B^j^ 
as a sum of two terms: 

Bijkix+) = - (a+7y(x+)) - - (5+7jfc(x+) - a+7ifc(x+)j . (106) 

The term containing a derivative of the delta function in the variation of the n-point function 
together with the first term in Bij^ give: 

ilmnixKix,) . . . C(X.)) ^ / (01^ • • • 1^(2-)' (- E ka^dk^^M^ 

Q Z 
s 

The form of this contribution is the same as for dilations. When they sum up the final result will 
be the same as for dilations but with a modified delta function 6{Y^ ki + q). This means that in all 
computations the n-point function has to be computed slightly off-shell. The two additional terms 
that contain a derivative of the delta function will give a total contribution proportional to 

~ dpJ{P)e%{q)q, J^ [k^.dk^^ - kajd^^^) M . (108) 

a 

This expression is equal to zero. In order to prove that, we can rewrite it taking into account that 
the amplitude is in general a function of scalar products ka ■ kb 

Y^ {kaidk,, - kajdk,,) M= Y ^7? ^ (^bjkai - kajktri) = . (109) 

a a,b=l "^(.^a • kb) 

Finally, the last part of the variation is 

^ . . . ^(27rf6{P)YiB^M^^)ka^dk^^dk^,Me^''-'^^-^^'--^- . 

(110) 
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After averaging over the long mode it gives the following contribution 

s 

Summing up all the nonzero terms and going to momentum space on the lefthand side of the 
equation, we get the generalized consistency relation for a background tensor mode: 

(7|Cfei • • • CfcJU = -\pii^) E 4(9i (^- + ^-(9"'- 4) - \{q- ka)dA dkjCk, ■ ■ ■ ky ■ (112) 

References 

[1] S. R. Bclibaliani, A. Dymarsky, M. Mirbabayi and L. Senatore, "(Small) Resonant non-Gaussianities: 
Signatures of a Discrete Shift Symmetry in the Effective Field Theory of Inflation," arXiv:1111.3373 

[hep-th] . 

[2] P. Creminelli, A. Nicolis and E. Trincherini, "Galilean Genesis: An Alternative to inflation," JCAP 
1011, 021 (2010) |arXiv:1007.0"027l [hep-th]]. 

[3] J. M. Maldacena and G. L. Pimentel, "On graviton non-Gaussianities during inflation," JHEP 1109, 
045 (2011) a rXiv: 1104.2846 [hep-th]]. 

[4] I. Antoniadis, P. O. Mazur and E. Mottola, "Conformal Invariance, Dark Energy, and CMB Non- 
Gaussianitv." larXiv:1103.4164 [gr-qc]. 

[5] P. Creminelli, "Conformal invariance of scalar perturbations in inflation," Phys. Rev. D 85, 041302 
(2012) ; arXiv:1108.0874l [hep-th]]. 

[6] A. Nicolis, R. Rattazzi and E. Trincherini, "The Galileon as a local modification of gravity," Phys. Rev. 
D 79, 064036 (2009) arXiv:0811.2197 [hep-th]]. 

[7] C. Burrage, C. de Rham, D. Seery and A. J. ToUey, "Galileon inflation," JCAP 1101, 014 (2011) 
|arXiv:1009.2497 [hep-th]]. 

[8] P. Creminelli, G. D Amico, M. Musso, J. Norena and E. Trincherini, "Galilean symmetry in the effective 
theory of inflation: new shapes of non-Gaussianity," JCAP 1102, 006 (2011) [arXi v:101l!3004l [hep-th]]. 

[9] C. Dcffayct, G. Esposito-Farese and A. Vikman, "Covariant Galileon," Phys. Rev. D 79, 084003 (2009) 
arXiv:0 901.1314 [hep-th]]. 

[10] J. M. Maldacena, "Non-Gaussian features of primordial fluctuations in single field inflationary models," 
JHEP 0305, 013 (2003) [arXiv:astro-ph/0210603] . 

[11] P. Creminelli and M. Zaldarriaga, "Single held consistency relation for the 3-point function," JCAP 
0410, 006 (2004) [astro-ph/04 07059|. 

[12] C. Cheung, A. L. Fitzpatrick, J. Kaplan and L. Senatore, "On the consistency relation of the 3-point 
function in single field inflation," JCAP 0802, 021 (2008) arXiv:0709.0295 [hep-th]]. 

[13] P. Creminelli, G. DAmico, M. Musso and J. Norena, "The (not so) squeezed limit of the primordial 
3-point function," JCAP 1111, 038 (2011) [arXiv:1106. 14621 [astro-ph.CO]]. 

[14] S. Weinberg, "Adiabatic modes in cosmology," Phys. Rev. D67, 123504 (2003). |astro-ph/0302326| . 

[15] X. Chen, B. Hu, M. -x. Huang, G. Shiu and Y. Wang, "Large Primordial Trispcctra in General Single 
Field hiflation," JCAP 0908, 008 (2009) tarXrv :0905.3494i [astro-ph.CO]]. 

25 



[16] F. Arroja, S. Mizuno, K. Koyama and T. Tanaka, "On the full trispectrum in single field DBI-infiation," 
Phys. Rev. D 80, 043527 (2009) |arXiv:0905.364n [hep-th]]. 

[17] C. Cheung, P. Creminelh, A. L. Fitzpatrick, J. Kaplan and L. Senatore, "The Effective Field Theory of 
Inflation," JHEP 0803, 014 (2008) arXiv:0709.0293 [hep-th]]. 

[18] L. Senatore and M. Zaldarriaga, "A Naturally Large Four-Point Function in Single Field Inflation," 
JCAP 1101, 003 (2011) [arXiv:1004.1201 [hep-th]]. 

[19] D. Seery, J. E. Lidsey, M. S. Sloth, "The inflationary trispectrum," JCAP 0701, 027 (2007). 
[astro-ph/0610210l . 

[20] D. Seery, M. S. Sloth, F. Vernizzi, "Inflationary trispectrum from graviton exchange," JCAP 0903, 018 
(2009). [ arXiv:0811.3934. [astro-ph]]. 

[21] X. Chen, R. Easther and E. A. Lim, "Generation and Characterization of Large Non-Gaussianities in 
Single Field Inflation," JCAP 0804, 010 (2008) jarXiv:0801.3295 [astro-ph]]. 

[22] L. McAllister, E. Silverstein and A. Westphal, "Gravity Waves and Linear Inflation from Axion Mon- 
odromy," Phys. Rev. D 82, 046003 (2010) |arXiv:0808.0706l [hep-th]]. 

[23] L. Leblond and E. Pajer, "Resonant Trispectrum and a Dozen More Primordial N-point functions," 
JCAP 1101, 035 (2011) [arXiv:1010.4565 [hep-th]]. 

[24] E. Pajer and M. Zaldarriaga, "A New Window on Primordial non-Gaussianity," larXiv:1201.5375 [astro- 
ph.CO]. 

[25] R. Scoccimarro, L. Hui, M. Manera and K. C. Chan, "Large-scale Bias and Eflficient Generation of Initial 
Conditions for Non-Local Primordial Non-Gaussianity," larXiv:1108.5512l [astro-ph. CO]. 

[26] C. Armendariz-Picon and E. A. Lim, "Scale Invariance without Inflation?," JCAP 0312, 002 (2003) 
arXiv:astro-ph/0307101 . 

[27] J. Khoury and F. Piazza, "Rapidly- Varying Speed of Sound, Scale Invariance and Non-Gaussian Signa- 
tures," JCAP 0907 (2009) 026 jarXiv:0811 .3633 [hep-th]]. 

[28] J. Khoury, P. J. Steinhardt, "Adiabatic Ekpyrosis: Scale-Invariant Curvature Perturbations from a Single 
Scalar Field in a Contracting Universe," Phys. Rev. Lett. 104 (2010) 091301. [arXi v:0910.2230l [hep-th]]. 

[29] J. Khoury and P. J. Steinhardt, "Generating Scale- Invariant Perturbations from Rapidly- Evolving Equa- 
tion of State," Phys. Rev. D 83, 123502 (2011) [ arXiv:1101.3548 [hep-th]]. 

[30] V. A. Rubakov, "Harrison-Zeldovich spectrum from conformal invariance," JCAP 0909, 030 (2009) 
larXivL0906.3693 [hep-th]]. 

[31] K. Hinterbichler and J. Khoury, "The Pseudo-Conformal Universe: Scale Invariance from Spontaneous 
Breaking of Conformal Symmetry," larXiv:1106. 14281 [hep-th]. 

[32] K. Hinterbichler, A. Joyce and J. Khoury, "Non-linear Realizations of Conformal Symmetry and Effective 
Field Theory for the Pseudo-Conformal Universe," arXiv:1202.6056 [hep-th]. 

[33] K. Hinterbichler, L. Hui and J. Khoury, "Conformal symmetries of adiabatic perturbations in cosmology" , 
to appear. 

[34] K. Hinterbichler, L. Hui and J. Khoury, work in progress. 

[35] W. Goldberger, L. Hui and A. Nicolis, work in progress. 

[36] L. Boubekeur, P. Creminelli, J. Norena, F. Vernizzi, "Action approach to cosmological perturbations: 
the 2nd order metric in matter dominance," JCAP 0808, 028 (2008). [arXiv:0806T0T6^ [astro-ph]]. 



26 



